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Abstract 

We study the non-continuous correction in the dynamics of drop impact on a solid 
substrate. Close to impact, a thin film of gas is formed beneath the drop so that 
the local Knudsen number is of order one. We consider the first correction to the 
dynamics which consists of allowing slip of the gas along the substrate and the 
interface. We focus on the singular dynamics of entrapment that can be seen when 
surface tension and liquid viscosity can be neglected. There we show that different 
dynamical regimes are present that tend to lower the singularity strength. We finally 
suggest how these effects might be connected to the influence of the gas pressure 
in the impact dynamics observed in recent experiments. To cite this article: L. 
Duchemin and C. Josserand, C R. Mecanique 333 (2012). 
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1 Introduction 



Drop impact is crucial in many multiphase flows ranging from raindrops to 
combustion chambers or ink-jet printing and it has become an emblematic 
problem of surface flows [1,2]. Depending on the impact parameters (drop 
diameter, velocity), fluid properties (viscosity, density, surface tension) and 
impacted surface (liquid deep or thin film, solid substrate), it can lead to 
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many different outcomes: spreading, rebound, prompt splash, crown splash, 
cavity and jet formation to cite the most famous ones [3]. Often, the influ- 
ence of the surrounding gas is neglected in the analysis because of the high 
density and viscosity ratios between the gas and drop liquid. Indeed, beside 
the entrapment of a gas bubble at impact due to lubrication effect underneath 
the drop [4,5] and some aerodynamic corrections to the corolla dynamics, no 
significant effects of the gas was noticed so far. However, the situation has 
suddenly changed recently in a striking experiments on drop impacts on a 
smooth solid substrate [6]: there, by changing only the operating pressure, 
they observe that splashes were suppressed as the pressure was lower below a 
critical level, emphasizing thus the crucial role played by the gas in the splash- 
ing dynamics. Although different scenari has been proposed to explain such 
effect, involving in particular gas compressibility [6], singular bubble entrap- 
ment dynamics [7,8], thin film skating [9] and film wetting dynamics [10,11], 
the surrounding gas influence remains yet a vibrant question of scientific de- 
bates. 

Of particular interest is the coupled dynamics between the drop and the gas 
underneath it just before the impact. In this case, it can be shown that the 
thin film air dynamics can be considered within the lubrication approxima- 
tion while the liquid viscosity can be neglected as far as thin liquid jets are 
not formed [7,9,12]. Then neglecting the surface tension one can see that a 
finite time singularity arises as a gas bubble is entrapped by the dynamics. 
This singularity behavior has to be regularized physically at least by the sur- 
face tension and the liquid viscosity but it has been argued that the resulting 
violent dynamics might be relevant in the splashing dynamics. Within the 
lubrication approximation for the gas film, a liquid jet is then formed that 
skates on the very thin (but non zero) gas layer. Eventually it has been shown 
experimentally that the liquid wets the solid substrate [10,11], something that 
cannot be explained within such lubrication approximation when surface ten- 
sion is present. In fact, different effects can be proposed to explain the liquid 
contact with the substrate when the gas layer is very thin: notably surface 
(Van der Waals for instance) interaction between the liquid interface and the 
substrate, interface fluctuations and/or surface roughness, and finite size (or 
rarefied gas) effect in the gas layer. In this paper, we focus on this latter case, 
that is when the gas layer thickness becomes of the order of the mean free 
path of the gas, leading firstly to a corrected lubrication equation for the thin 
film. In the next section, we recall the general scaling argument obtained for 
drop impact on a solid substrate using the classical lubrication equation for 
incompressible fluids. Then in section 3 we introduce the correction when the 
gas thickness is of the order of the mean-free path. Finally we discuss in sec- 
tion 5 the properties of the singularity in this case, in the absence of surface 
tension and we draw some perspectives for this work. 
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Figure 1. Sketch of the impacting drop and successive profiles for 
We = piDU 2 h = 95 and St = 1.35 x 1(T 3 . 



2 Scaling analysis 



We consider the impact of a liquid drop of diameter D with a vertical velocity 
U on a smooth solid substrate (Cf. Fig. 1). The liquid and surrounding gas 
densities are noted pi and p g , their dynamical viscosity m and \i g respectively, 
the surface tension 7. Considering the high Reynolds and Froude numbers of 
the drop for typical experimental conditions 

Re=^ Bo U 



pi VgD 

we assume that we can neglect the liquid viscosity and the gravity in the 
dynamics that is driven by the cushioning of the air film beneath the drop at 
impact. 

Therefore, considering incompressible fluids we obtain the following set of 
equations describing the evolution of the drop surface h(r, t) in axisymmetric 
geometry when it is approaching the substrate: 



(dn) d t ip + \ Vip 2 -\- — h — n = c(t), (i) 
2 pi pi 

(dn) d t h = -!—d r (rh 3 d rP ), (2) 
iZrp,g 

(dn) d t h = d z (f — d r (fd r h, (3) 

(n) Aip = 0, (4) 

where ip is the liquid velocity potential (the velocity in the liquid is u = V<p), 

and n is the drop volume, dn its interface. The first equation (1) is the 
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Bernoulli equation valid in the liquid and written at the interface, while the 
last one (4) is the incompressible condition for such potential flow. The inter- 
face motion is described by the advection equation (3) while the lubrication 
equation (2) allows the determination of the pressure p(r, t) in the gas film. 
The 1/12 numerical prefactor in this equation was obtained by considering no 
slip condition for the gas velocity both on the solid substrate and on the fluid 
interface. This system of equations can be solved numerically using a bound- 
ary integral method and we resort with numerical simulation on the interface 
only [9] . Notice that lubrication is only valid where the slope of the interface 
is small enough and additional terms should be considered otherwise. 
As the drop is approaching the solid substrate, cushioning of the gas leads to 
high pressure gradient beneath the drop. The high pressure created at the bot- 
tom of the drop deforms its shape so that a gas pocket is entrapped around the 
impact center. Simple scaling arguments can help to estimate the typical sizes 
of this entrapped gas bubble. Indeed, considering the time t = as the time 
of impact in the absence of air, one can estimate the typical vertical H and 
horizontal R scales of the drop deformation as H(t) ~ Ut and R(t) ~ \JDH 
thanks to a geometrical argument based on the intersection between the falling 
drop and the substrate. Introducing these scalings in the lubrication equation 
yields the following scaling for the lubrication pressure in the gas layer: 

On the other hand, considering the impact pressure in the liquid needed to 
deviate horizontally a volume of liquid of typical size R (this argument was 
first given in [13]), one obtains : 

Thus, the lubrication pressure Pi is strong enough to deviate the liquid until a 
critical height H* where the liquid somehow has to touch the solid substrate, 
that is when Pj = Pi, yielding: 

/ \ 2/3 
H*~D =PSt 2 / 3 

introducing the Stokes number as the inverse of a Reynolds number balancing 
the liquid inertia with the gas viscous effects: 

St- ^ 

Therefore, one expects the liquid to contact the substrate for t* ~ St 2 / 3 D/U 
entrapping a gas bubble of height H* and radius R* ~ DSt 1 ^ 3 . Notice that 
in practical (experimental) situations, the Stokes number is very small: for 
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instance for a 2 mm diameter droplet of water impacting at a velocity of 1 
meter per second, we obtain St ~ 1CT 8 , so that the entrapped bubble radius 
is only a few thousands of the drop diameter. 



However, it has been observed in numerical simulations of the set of equa- 
tions (1,2,3,4) that this contact would arise as a finite time singularity in the 
absence of surface tension (taking 7 = in equation 1) [12,7,9]. Such singu- 
larity exhibits a divergence of the pressure following P ~ h^J^ and of the 
interface curvature k ~ h^ in which can be explained following [9]. Writing 
the set of governing equations in the frame moving radially with the geomet- 
rical intersection between the falling sphere and the substrate and using the 
dimensionless variables defined by 

„ r - R(t) „ z Ut ~ h(r,t) _ p _ cp 

x = — , z — — , t — — , n{x,t) = — — — , p = — — , and ip 



D ' D' D ' v ' ' D ' ' p^ 2 ' Y UD' 

the following system of equations is obtained: 



d~ t Cp - kdiCp + U/f +P = C(t), (5) 

dfh - Rdji = 1 ((x + St 1/3 )h 3 dxp) , (6) 

* 12St(x + St 1/3 ) ^ ' FJ ' y ; 

dfh — Rdxh = d z (p — d r (pd r h, (7) 
A(p = 0. (8) 

The dimensionless radius R = \J H / D gives for the singularity a dimensionless 



radial velocity R ~ yD/H. 

Thus we develop the former set of equations near the singularity where R ~ 
St -1 / 3 (we drop the ~ thereafter for the sake of simplicity) and assume that 
the time derivatives are subdominant compared with the Rd x terms, in good 
agreement with the numerics [9]. Then seeking for a self-similar structure for 
the interface near the singularity in the form: 

h(x,t) = h min f(—), p(x,t) = P (t)g( — ) and (p(x,t) = (p (t)$(j^, ^y) 

we obtain the following system of equations at the interface at the dominant 
order in the self similar variable £ — x/l(t) and x — z/l(t): 

- St-/ 3 ^^ + (V$) 2 + P (t)g = C(t) (9) 

bt i(t) ! ~ i2st i{ty u 9 > 
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where the prime stands for the derivative of the function to the variable £. 
The lubrication equation (10) gives the following relation: 

and from the interface dynamics eq. (11) one can see that two situations have 
to be considered: 

- h m in(t) <C l(t) leading to the observed numerical scalings since <po(t) ~ 
St~ 1/3 h min (t) that gives P (t) ~ St' 2/3 h min (t) /l(t) and thus: 

l(t) ~ St- 2/3 h mm (tf 2 , P (t) ~ h^t)- 1 / 2 and k ~ St 4 / 3 /i min (*)- 2 . 

and this regime is found to be valid for thick films, h min ^> St 4 ^ 3 . 

- h min (t) 3> l(t), which gives (po(t) ~ St _1 / 3 Z(i) and yielding: 

P (t) ~ St" 2 / 3 , Z(f) ~ St- 4 / 3 /i mm (t) 2 and k ~ St 8 / 3 /i m in(0 -3 - 

Similarly, this regime is valid for thin films, h m i n St 4 ^ 3 . Finally, notice 
that this second regime has never been reached numerically so far because 
of the small values of the Stokes number considered. 



3 Influence on the mean free path 



It is interesting to observe that under typical experimental conditions the 
Stokes number is very small so that the system of equations used has to be 
questionned. In particular, when the typical air layer becomes of the order of 
the mean free path A one has to consider rarefied gas correction to the gas 
dynamics. This can be quantified by the Knudsen number Kn defined as the 
ratio between the mean free path and the typical air layer thickness which 
gives for the entrapped bubble described above: 

and one expects rarefied gas effects to enter into account for Kn ^> 0.01. From 
kinetic theory, we have that the mean free path is related to the gas pressure 
P g through: 

A = k B T 

^2lXd 2 Pg' 
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where T is the temperature, ks the Boltzmann constant and d the typical size 
of the atoms of the gas. For ambient temperature T = 300 K and ambient 
pressure P g0 = 10 5 Pa the typical mean free path in the air is A = 70 nm and 
one can write the simple relation: 

A = PgQ 

Ao Pg 

In the situation of incompressible fluids (gas and liquid) considered here, A 
is a constant that is only a function of the ambient pressure and it does not 
formally depends on the dynamical pressure used in the equations. In fact, 
such approximation is valid as long as one can neglect the gas compressibility 
which should be accounted for otherwise. The effect of the gas compressibility 
has been studied in details in [7,8,14] using a low Mach number lubrication 
approximation but no crystal clear mechanism involving directly the gas com- 
pressibility could be identified in the drop splashing, beside some complicated 
dependance of the dynamics. Therefore, in order to disentangle the influence 
of the rarefied gas correction from the compressible influence, we keep the 
incompressible limit in the dynamical equation, considering only the influence 
of the Knudsen number in the boundary conditions as explained below. Fi- 
nally, as the rarefied gas situation is intrinsically for compressible fluids, an 
extension of this work might have to consider compressible effects. 

Since A, and thus Kn, increases as the gas pressure decreases, one can expect 
that the correction due to the mean free path of the gas has to be considered 
when the gas pressure is lowered. Indeed, taking the experimental conditions 
of [6] where the striking influence of the air pressure was first illustrated, one 
finds for atmospheric pressure Kn = 0.66 so that correction due to the rarefied 
gas configuration has to be investigated. 

For such relatively high Knudsen numbers, a simple way to account for the 
correction due to the rarefied gas context is to introduce a slipping velocity for 
the gas so that the no-slip boundary condition at a solid interface transforms 
into the Navier-slip condition: 

,du t 
on 

where n is the normal direction at the interface and u t the tangential velocity. 
This condition comes from the fact that at the level of the mean free path the 
no-slip condition is meaningless and cannot be imposed. Solving the Stokes 
equation between the solid substrate and the drop interface located at z — 
h(r,t) and assuming that the Navier-slip condition applies on both sides we 
obtain the following relation for the radial velocity u in the gas layer, under 
the thin film approximation where the interface slope is supposed to be small: 

u(z,t) = -^—d r p(hz-z 2 + \h) , (12) 
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so that the lubrication equation (2) becomes: 



d t h = 



— (d r {rh 3 d rP ) + 6Xd r {rh 2 d r p)) . 



(13) 
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Therefore, for film thickness /i< A one expects that the second term in the 
lubrication dominates so that the singularity features should be changed. 



4 Scaling in the rarefied gas limit 

The general dynamics of air cushioning during drop impact will thus be gov- 
erned by equation (13) and different dynamical regimes will dominate depend- 
ing on the air layer thickness. When the air layer is always much thicker than 
the mean free path A then the former equation (2) will be valid and the scaling 
obtained in section 2 will be observed. On the other hand, when /i < A the 
lubrication dynamics will be dominated by the new term introduced in (13) 
so that the following equation should be investigated: 



where the usual lubrication term is neglected. In between and when the air 
layer is varying within these limits, one has to study the general equation 
(13). Moreover, two different situations can be identified: firstly, Kn C 1 so 
that the bubble entrapment still arises around the same radius DSt 1 ^ 3 than 
above, although the dynamics of the film cushioning at the singularity neck will 
eventually be governed by the slipping equation (14). On the other hand, for 
Kn > 1, already the dynamics of the bubble entrapment has to be considered 
within the rarefied gas limit of eq. (14) so that the radius of the bubble itself 
(and thus the radial velocity of the interface near the singularity) is changed. 

4-1 Low Knudsen numbers case 

First of all, we thus consider that the gas bubble entrapped has the same 
features than before. This assumption is reasonable for small enough Knud- 
sen number and the bubble formation is dominated by the usual lubrication 
equation. However, the singularity dynamics itself has to be studied within 
the rarefied gas limit since the local air layer becomes much thinner than the 
mean free path. Then we can write the dynamics in the frame moving with the 
former singularity velocity USt^ 1 ^ 3 . Using the same dimensionless units and 




(14) 
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then seeking for similarity solution using the same change of variables than 
before, we obtain at the dominant order near the singularity: 



- St-" 3 ^|a £ * + (V#) 2 + P„(t) 9 = C(t) (15) 

-^'^r - ^MS^ 1 (/v + «^/v)' (i6) 

^t^SKW (17) 



Notice that A is now the mean free path made dimensionless using the drop 
diameter. Now, two regimes can also be identified due to the two terms in the 
right hand side of the lubrication equation (16). 

- h min {t) ^> A, the pressure relation due to the dominant term in the lubri- 
cation equation remains: 

p (t)~st 2 / 3 l[t) 



hminit) <C A, the pressure relation is determined by the other term in the 
lubrication equation, yielding 

P (t)~St 2 / 3 Kt) 



A/imin(^) 



Finally, we resort with two dynamical scenari for the singularity dynamics, 
depending on the ratio between the mean free path A and the critical thickness 
St 4 / 3 separating the two regimes of the singularity without slip condition, 
namely: 

- A > St 4 / 3 , two dynamical regimes follow. First, when h min 3> A, then the 
"usual" scaling are valid: 

lit) ~ St- 2 / 3 h min (t) 3 / 2 , P (t) ~ h min (t)- 1/2 and k ~ St^ 3 h mm (t)- 2 - 

then it is followed by another regime as h min decreases. When h min <C A, 
the following scalings are obtained: 

1 St 4/3 
l(t) ~ y/\h min (t)St~ 2/3 , P (t) ~ —= and k 



and we remain thereafter within the configuration where h m i n (t)/l(t) ~ 
St 2/3 /V\ < 1. In this case, the regime where h min (t) ^> l(t) is not present. 



9 



A < St 4/3 , then for h min > St 4/3 > A the "usual" scaling holds: 

lit) ~ St- 2/3 h min itf 2 , P it) ~ /i mm (t)~ 1/2 and « ~ St 4/3 /i min (t)- 2 . 

It is followed by a regime where St 4 ^ 3 ^> h min ^> A which gives the second 
scaling obtained in the beginning 

Poit) ~ St' 2 / 3 , Z(f) ~ St- 4/3 /i mm (t) 2 and k ~ St 8/3 /i mm (t)- 3 . 

Such regime is also followed by a new regime when St 4 / 3 ^> A ^> h min 
yielding: 

Q+8/3 

lit) ~ A/w(t)St~ 4/3 P (0 ~ St~ 2/3 and k 



A 2 ^l m jn(^) 

^.2 Full rarefied gas limit 



In the experiments described in [6], the Knudsen number based on the typi- 
cal bubble thickness obtained with the usual lubrication equation is of order 
one and it increases when the pressure is lowered. Therefore, the rarefied gas 
effect has to be accounted for in the dynamics of the bubble entrapment and 
the above scalings fail already for the size of the entrapped bubble. In this 
regime, we have A 3> -DSt 2 ^ 3 3> h(r, t) everywhere and the second term in the 
lubrication equation - proportional to A - is dominant. Therefore, the scaling 
for the lubrication pressure reads : 

y g R 2 U 
' XH 2 

Balancing the impact pressure with this lubrication pressure, one obtains for 
the typical bubble height : 

It needs to be emphasized that we obtain here a bubble height H* and radius 

R* = v/M-DSt^, 
A 

that depend on the external gas pressure through A! Notably, the bubble 
radius is proportional to the external pressure and therefore decreases when 
this ambient pressure is lowered. Considering again that the singularity follows 
the geometrical radial scale R ~ y/DUt, we obtain that the radial velocity of 
the singularity yields: 
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Again, we can assume that the time derivatives are subdominant compared to 
the Rd x terms. Then, after writing the dynamical equation in dimensionless 
form and seeking for self-similar solutions, we obtain the following system of 
equation at the interface at the dominant order in the self similar variable 
f = x/l(t) and X = z/l(t): 



i¥ v+ w (w), ; ,M "' =c( " (18) 

A h min (t) , _ 1 h min (t) AP (t) / 2 ,y 

St Z(t) 7 ~ 2St /(t) 2 ^ 5 ) ^> 



st z(*) J i(t) V 

where we have assumed that A ^> h is true everywhere in the bubble region 
so that only the slipping term in equation (19) is present. 

Therefore, balancing the two terms in the lubrication equation, we obtain for 
the pressure : 

p m 

Po{t) h~TtY 

Then, we get again a priori two different regimes, depending on the value of 
hmin compared to the value of l(t) : 

- h min (t) <l(t), leading to (p (t) ~ Xh min (t)/St that gives : 

l(t) ~ \h min (t)St-\ P (t) ~ ASt- 1 and k ~ St 2 A- 2 /i m , n (0 _1 , 

and this regime is found to be valid when A 3> St, which is always the case 
since this regime was found to be valid for A ^> St 2 / 3 ^> St (since St 1 in 
the experiments). 

- h min (t) 3> l{t), would give po(t) ~ A/(t)/St and yields: 

lit) ~ \ 2 h min (t)St~ 2 , P (t) ~ A 2 St- 2 and k ~ St A \- 4 h mm (t)-\ 

As explained above, this regime cannot be observed in realistic experiments, 
since it implies A St, and therefore that the bubble height is much bigger 
that the radius of the drop. 
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5 Discussion 



We have obtained the correction due to the account of the rarefied gas limit in 
the singularity dynamics leading to a bubble entrapment during drop impact. 
In particular, we have shown that the pertinent regime for the experiments, 
where the splashing properties are influenced by the gas pressure, corresponds 
to a situation where the bubble size reduces with the gas pressure. The ques- 
tion of vibrant scientific debates is whether and how one can relate this singu- 
larity to the splashing dynamics of drop impact on solid substrate. Indeed, it 
has been often argued that such violent singular behavior is at the heart of the 
onset of splashing so that the regularization of the singularity would be crucial 
to determine the splashing properties [7,8,9,14]. In particular, surface tension 
on one side and liquid viscosity on the other side will both regularize the sin- 
gularity, leading to a rapid jet propagating along the solid wall [9]. However, 
such effects could not explain the experimental results of [6] if only the usual 
(without the rarefied gas correction) incompressible lubrication dynamics was 
considered for the air layer. This is why compressible effects have been often 
investigated, in the framework of the low Mach number compressible lubrica- 
tion equation [7,8]. Finally, while within the lubrication approximation, the jet 
skates on the thin gas layer, it has to be noticed that experimentally the liquid 
eventually wets the solid substrate [11,10] through a dynamical rupture of the 
gas layer that still needs to be elucidated. Here we would like to emphasize 
that the rarefied gas effects accounted in the incompressible limit offers an 
alternative physical mechanism to explain the air pressure effect in the drop 
impact, as suggested also in the compressible model [14]. In particular, we 
have demonstrated here that the entrapped bubble size depends strongly on 
the external gas pressure, a result that has to be investigated in future works 
since it changes the jet features and consequently the air drainage dynamics. 
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